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Solutions are obtained for several axially symmetric plane strain problems involving a hollow 
circular viscoelastic cylinder. The cylinder is assumed to be subject to planar, axially symmetric body 
force and temperature fields. Displacement, traction, and mixed boundary conditions are considered. 
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1. Introduction 

In this paper we study a hollow thick-walled right circular cylinder composed of an isotropic 
homogeneous viscoelastic material whose response to temperature is thermorheologically simple. 

We solve the axially symmetric plane strain problem for this cylinder subject to three types 
of boundary conditions: 

(i) zero displacement at the inner and outer walls, 
(ii) zero traction on both walls, and 

(iii) zero traction inside, zero displacement outside. 

The latter mixed problem was chosen to suggest the solid fuel core of a rocket. Our analysis 
can be easily modified to deal with the case of a solid cylinder or with nonvanishing boundary 
conditions. 

Many problems in isothermal viscoelasticity are readily solved by means of Laplace trans- 
forms. This technique has also been applied with some success in thermoviscoelasticity, for 
example to sphere problems by Muki and Sternberg [l] 1 and Sternberg and Gurtin [2], and to 
a contact problem by Graham [3]. 

In fact Valanis and Lianis [9] use Laplace transforms to obtain the solution of a thermo- 
viscoelastic cylinder problem under the assumption, unnecessary in the present paper, of a constant 
hereditary Poisson ratio. 

Our method is fundamentally different from the above. It is a generalization of that used in [4] 
to obtain existence for the displacement problem of isothermal viscoelasticity. The analysis in [4] 
led to a Neumann expansion representation for the viscoelastic displacement field which involved 
the inverse E~ x of an elastic operator E. Also, convergence of the expansion depended on a Schauder 
estimate for E. In order to assure existence of E~ x and of the estimate we were obliged in that 
paper to confine ourselves to the displacement problem. For the solution of the cylinder problem 
we obtain an analogous Neumann expansion. However, due to the extremely simple geometry 
of the present problem an explicit representation of E~ x is easily derived for a variety of boundary 
conditions. Thus, our solution is altogether explicit, and moreover, estimates which suffice for 
the convergence of the expansion are immediate. 

Finally we remark that no attempt has been made to expound on the derivation of the con- 
stitutive equations set forth below. Thorough expositions of the theory of linear viscoelasticity and 
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thermorheologically simple materials are given in the papers [5, 6] of Gurtin and Sternberg. If 
one accepts our field equations, then our analysis can be followed without consultation of any 
of the above references. 

We plan a sequel to this paper in which the effects of ablation are considered and numerical 
examples are given. 

2. Statement of the Problems 

Consider a hollow circular cylinder with inner radius a and outer radius b. If Xi, jc 2 , and x 3 
are orthogonal Cartesian coordinates oriented in such a way that the # 3 -axis coincides with the 
axis of the cylinder, then a set of cylindrical coordinates p, <p, zis defined by 

Xi = pcos(p, X2 = ps'm<p, Xs = z 
where 

0^p<oo, 0^<p<2iT, — 00< Z <00. 

The cylinder is composed of isotropic homogeneous viscoelastic material subject to a tem- 
perature field 2 T(p, t), where t stands for time, and a body force vector 

F=(F P ,F,,F.) 
such that 

F z = F* = 0, F p = F(p,t). 

We assume that the displacement field 

u= (u p , Wp, u z ) 
is axially symmetric and planar, i.e., 

u z =u<p = 0, u p =u(p, t). (2.1) 

Let <j pp , (Tw> • • • and e pp , e^, . . . denote respectively the components of stress and strain 
in cylindrical coordinates. Then from (2.1) and the usual strain-displacement relations for these 
coordinates ([7], p. 183) it follows that 

6 P p(p, t)=— (p, 0, 6w(p, t)= " , 6^ = 

dp p 

V €pj> €<fz €zp 

Let G\{t) , Gz(t) be respectively the relaxation functions in shear and isotropic compression 
corresponding to the base temperature To. We assume that the cylinder has been at rest at the 
reference temperature To up to time £ = and that the material with which we are dealing is ther- 
morheologically simple. 

In order to express concisely the plane strain stress-strain relations for a thermorheologically 
simple materia] and to avoid needlessly detailed equations in the sequel, we introduce the following 
operator notation. Given any differentiable function G(t) defined on [0,oo) we associate with it 
the thermoviscoelastic operator Vg by defining 3 

V G f(p, t) = G(Q)f(p, t) + £ GW ^ <p{T(p,T))dr)ip{T{p, f'))/(p, t')dt'. (2.3) 



2 In this paper, we have not attempted to determine the minimum sufficient smoothness of the given data. All given functions are assumed to be several times 
continuously differentiable. 
dC 



du 
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Here, Tip, t) is the temperature field introduced above, and <p(T), the so-called shift factor, is a 
prescribed function which governs the change in material properties induced by a temperature 
field. The function <p has the properties 

<p(T )= l,<p(7) >0,^>0. (2.4) 

Notice that for C(0) ¥" 0, V G is a Volterra-type integral operator of the second kind. We shall assume 
that all the thernioviscoelastic operators which we encounter are of the second kind. 
It is also convenient to define 

M{t) = h(G 2 (t)-G l {t)) (2.5) 

d(0)=y, M(0)=/x (2.6) 

g(pj,t')=G\^(y I ip{T{p,r))dr\ip[T(p,t f )] (2.7) 

m(p, t, t')=M {x) ( I <p(T(p, T))dj\ (p[T(p, t')\ (2.8) 

so that V(; x and Vm take the form 

V ( ;J(p, t)=yf{p,t)+ f'gip, t, t')f(p, t'W (2.9) 

Jo 



V M f(p, t)=pLf{p, t)+ m(p, t, t')f(p, t')dt' . (2.10) 

Finally we introduce the function 

<4>(p, t)=a {) V ( ,fi(p, t). (2.11) 

The function B is defined by 

B(p, t)=— \ a(T)dT\ a () =a(T {) ) 

«o J T 

and is known to thermoviscoelasticians as a "pseudotemperature." The prescribed material func- 
tion ol{T) is called the coefficient of linear thermal expansion. 

The thernioviscoelastic stress-strain relations appear in many papers, for example Sternberg 
and Gurtin [2]. Written in terms of Cartesian coordinates and with the present operator notation 
they take the form 

(Tijip, t)= V Gl €ij(p, t)+8ijVM€ kk (p, t)- 0;/I>(p, t). (2.12) 

In this equation, Latin subscripts take the range 1,2,3 and summation over repeated indices is 
implied. In the case of plane strain, Sternberg and Gurtin [2] apply Laplace transforms to solve 
for cr zz in terms of the other stresses in the form 



where 4 



ir/defined on [(),*) f*ir))= * e^ l f(t)dt. 

Jo 



crzz = V-zio-pp + o- w ) - ou,Vf3 (2.13) 

p* M= 3Gf(7j)G 2 *(r?) 
[V ' Gf(ri)+2G$(j)) 

- ni|) . -cm-em , 

W 7,(C I *(7,) + 2G 2 *(TJ) 
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From equations (2.12), which are invariant under a change of coordinate systems, and (2.2) 
we obtain the plane strain axisymmetric stress-strain relations 



OVv 



-^■(a7) +V -(g + |)-*- 

=M-:K(i+;)-* 



(2.14) 



(Tp<p — (7<p z — <J zp — 0. 

Due to (2.13 and (2.14), the quasistatic stress equations of equilibrium ([7], page 184) reduce to the 
single equation 



d 1 

— o-pp-\ — (op P — aw) + F = 0. 

dp p 



(2.15) 



If we substitute from (2.14) into (2.15) and apply the notation (2.9), (2.10) then the quasistatic 
displacement equation of motion takes the form 



^(p, {) ^^^ 



dp 



dp 



( u(p, t) \ 



■— p 

y+M Jo 



(g(p, t, t')+m{p,t, t')) 



>^im 



.(§£. 



+ {¥p^ t ' tn T P {p ' t ' t)) Tp^ t) 



dm u(p, t) 

+ T~ (P, t, t ) 

dp p 



\dt' = H(p,t) (2.16) 



where 



^•*>-^((0)^')-^4 



Three types of boundary value problems are considered, namely, 
(i) The displacement problem (first boundary value problem) 



u(a, t)=u(b, t)=0At>0) 



(ii) the traction problem 



(2.17) 



(2.18) 



(Tp P (a, t)=(Tpp(b, t)=0, and 



(2.19) 



(hi) the mixed problem 



0=cr pp (a, t)=u{b, t). 



(2.20) 



Our task then is to find a function w(p, t) which satisfies (2.16) for (p, t) in fa, b] X [0, oo) subject 
consecutively to each of the boundary conditions (2.18), (2.19), (2.20) where cr pp is defined in terms 
of u by (2.14). Stresses are then computed via (2.13), (2.14). 
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3. The solution 

We define the differential operator £ and the integrodifferential operator B as follows: 

Bfip, t)=-J— \ {(g(p, t, t')+m(p, t, t'))Ef(p, t') 
T + M Jo 



dp P J 



Then (2.16) reduces to 



Eu + Bu = H. (3.3) 

Our plan of attack is as follows: consider the model equation 

Ef+Bf=P (3.4) 

for which we seek a solution /(p, t) in fa, b] X [0, ») subject to the boundary conditions 

Mafia, t) = M h f(b, t)= 0(t> 0). (3.5) 

Here M (l and Mi, are some sort of linear integrodifferential boundary operators. Let E~ l be that 
operator which assigns to a given function t//(p, t) a unique function/ which satisfies 

£/=* (3.6) 

in [a, 6[ X [0, oo) and (3.5) for t ^ 0. Then clearly we obtain a solution/ for (3.4), (3.5) by first finding 
a solution i// for the equation 

ifj + BE~^ = P (3.7) 

and then setting 

/=£-'</<. (3.8) 

The advantage to be gained by dealing with (3.7) rather than (3.4) lies in the fact that (3.7) is a 
sort of generalized integral equation and may be solved by the method of successive approximations. 
This decomposition of a viscoelastic problem into an elastic problem and an integral equation was 
introduced by Fichera in [8]. 

The formal solution of (3.7) is 

ifj=(I + BE~ 1 ) 1 P 

i//=J (-IHBE-^yP (3.9) 

so that 

f= J (-iyE-HBE-iyP. (3.10) 

n = 
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We shall now set forth sufficient conditions for the uniform convergence of the series (3.9), 
(3.10) in any rectangle [a, b] X [0, t] (t finite). 
To this end, we define the norms 5 



(t) = l.u.b.{\Hp, t)| : (p, r)e[«, b] X [0, t]} 



i(t) 



(«) + 



(0 



(3.11) 
(3.12) 



and postulate the existence of a nonnegative nondecreasing function N(t) such that for any smooth 
function 0(p, t) and t^O, 



\\E-Wl(t)^N(t)\\0\\(t). 
It is clear from (3.2) that for smooth g and m, B satisfies an inequality of the form 

\\Bh\\(t)^Q(t) j\\\EhUt') + \\k\\i(t'))dt'. 

For example, one may take 

Q(t)= ,max \l.u.b.{\gfr,T,s) +m(p, r, s)\ : (p,t,s) e[a, b] X [0, t] X [0, r]}, 

|y+Ml I 



(3.13) 



(3.14) 



l.u.b. 






+ 



5m , v 



: (p, r, 5)e[a, 6] X [0, j] X [0, t] 



The function ^ so defined is clearly nonnegative and nondecreasing. 

By replacing h by E~ l 6 in (3.14) and applying (3.13) we obtain the estimate 



\\BE-^e\\(t)^Q(t)(i+N{t)) £\\e\\(t')dt'. 



(3.15) 



From (3.15), it follows by a standard iteration argument that 



\\(BE-*)*e\\(t) 



(Q(t)(l + N(t))t)» 



1(0. 



(3.16) 



Thus the series (3.9) converges uniformly in any rectangle [a, b] X [0, t]. Since (3.13) implies that 
E' 1 is continuous with respect to the uniform norm in such a rectangle, it follows that (3.10) also 
converges uniformly in [a, b] X [0, t] and that interchange of E~ l and the sum sign is permissible. 
Thus, sufficient conditions for the uniform convergence of (3.9) and (3.10) are the existence of the 
nondecreasing bounding functions /V(0 and Q(t). The boundedness of Q follows if we assume that 
d, G 2 , (f, and T are sufficiently smooth. The function N(t) depends on the latter functions together 
with the boundary conditions. 

Each of the boundary value problems (i), (ii), (iii) may be reduced to a problem of the form 
(3.4), (3.5). However, the operators M (l and M^ appearing in (3.5) will be different in each of the three 
cases. Thus, we must derive three different inverse operators J^j" 1 , Ef 1 >E$~ 1 corresponding respec- 
tively to (i), (ii), and (iii) and find nonnegative, nondecreasing functions Nu N2, N$ such that for any 
function <p and t ^ 0, £~j and Nj satisfy (3.13). Then by the above theory, the solution of each of the 
three problems will have a uniformly convergent representation of the form (3.10). 



s and h denote arbitrary smooth functions of p, t. 
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(i) The Displacement Problem 

Let E^ l 6 be that function h which satisfies 



Eh(p, t)=d(p, t) 



(3.17) 



in [a, 6]X[0, oo) 

and 

h(a, t) = h{b, t) = (3.18) 

for t ^ 0. By elementary techniques we obtain from (3.1), (3.17), and (3.18) the representation 

**» '>= \ />• » *?* *-#3 r e <* * *r* *>■ ™ 

It is obvious from (3.19) that a constant N\ depending only on a and b exists such that 

||/sr'€>ll.(') ^,11911(0. (3.20) 

Thus the displacement problem (2.18), (3.3) has the solution 

z/=jr (-\)"E-'(BE-')"H, ( 3.21) 

w.=0 

where Ef 1 is given explicitly by (3.19). 
(ii) The Traction Problem 
For this problem, it is convenient to modify the operator notation (2.9), (2.10) as follows: 



V<;,[/u/, t')](i)=yf(d, t)+Pg(c, /. t')f(d, t'W 



V' M \f(d, t')](t)=nf(d, ')+ rn(c, t, t')f(d, t')dt' ■ 



(3.22) 



/: 



Then by (3.3), (2.14), and (2.19) the traction problem takes the form 

Eu + Bu = H 



0(r/, t)=V Gi 



du , ,, 
— (a, t') 



a>(6,0=^ 



du 



(b, f) 



if)+V% 



(»+n 



du . ,. . u(a, t' 

— (a, t + 

dp a 



du ^ , u(b, t') 
^^t^—b— 



(t) 



w. 



We seek a solution u for this problem in the form 

u = v+w 



(3.3) 
(3.23) 

(3.24) 
(3.25) 



where w is to be chosen in such a way that v satisfies a problem of the form (3.4), (3.5). Plugging 
(3.25) into (3.3) we obtain 



Ev + Bv = H—Ew— Bw. 
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(3.26) 



Also we demand of w that it satisfy the boundary conditions 



<t>{a,t) = V« Gi 



*(b 9 t) = V Gi 

Then u denned by (3.25) will satisfy (3.23), (3.24) provided v satisfies 

dv 



dw , M 
— (a, O 

Lap J 


(*) + Vi 


dw t 

Lap 


t w(a, t') 
a 


(t) 


(3.27) 


aw 

— (6, t') 

Lap J 


(t)+n, 


dw ,. 

Lap 


w{b, t') 
} b . 


(*)• 


(3.28) 



0=F« Gi 



0=F G] 



dp 
aw 



(a,0 



(O + Fi 



ai; , ,v . v(a, t') 

— (a, t) +— 

dp a 



w 



»] 



(b,t')\ (t) + vi, 



— (*,«') + 

dp o 



6 J 



(«). 



The function w is by no means uniquely determined. We assume one of the form 

w(p, t)= {p-a)* <pi(t)+ (p-bV <pi(t) 



(3.29) 
(3.30) 

(3.31) 



where the functions <pi, <p2 are to be determined by a substitution of (3.31) into (3.27) and (3.28). 
In this way we arrive at the equations 



3>(a, = 



2(y+p)(a — b) + p 



(a-bV 



<pt(t) 



*(6,0 = 



+ I 1 2 (a — b)[g{a, t, t') + m{a, t, t')] + m(a, t, t') — — | <p 2 {t')dt' 

(*-o)«- 



2(y+/»)(6-o) + M 



<pi(0 



/:t 



+ ]2(6-a)[g(6, *, f')H-/n(6, t, ^')] + m(6, t, *') 



„ (b~ 



■«) 2 1 / »w» 



Since these are, in general, Volterra-type integral equations of the second kind their solution 
presents no difficulties. We may thus regard w as known, so that the traction problem has been 
reduced to that of finding a function v which satisfies (3.26), (3.29) and (3.30). To this end, we define 
E^ x to be that operator which assigns to a given function 0(p, t) the function h(p, t) which satisfies 
(3.17), (3.29), and (3.30). In order to derive E% *, we write the general solution of (3.17) in the form 



h(p, t) = - p f'diy, t) py^) dy+ Cl {t)p- 



Ci{t) 



Then 



dh 

dp 



(p 



.0— W' 

P 2 J a 



Q e{y,t)(£-Ady+\* 'e^Ddy+cAt)-^-- 

P' J a \ 2 / J a P~ 



(3.32) 



(3.33) 



The functions C\(t), Ci{t) are determined by a substitution of the expressions (3.32), (3.33) into 
(3.29), (3.30). They are governed by the equations 



0=Vg i+2U [c 1 (t')](t)-- 2 Vg l [c i (t')](t) 



r(0 = Pfc I+ ufM*')](*)-£rfc 1 [e*(«')](0 



(3.34) 
(3.35) 
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where 



r«)=n,[£ />,.') 



6 2 -y 2 



rfy 



w-n 



G x +M 



^ a 



6{y,t')dy 



(0. 



(3..%) 



Except in singular cases, which we choose to disregard, all of the operators of the form V G are 
Volterra-type integral operators of the second kind. Thus, they all have inverses, V~ l , which are 
operators of the same kind. Hence the system (3.34), (3.35) is readily solved, and we have 






7 2 ^ d+2M^ G, U2 V G t +2M r G, 



^WHOIM 



(3.37) 



c 2 (0 = 



— F->n J/" F" 1 * J 76 

a 2 ' Gt+2M y G, £2 ' G,+2M y G, 



PsiWT (»')](«). 



An explicit representation for E 2 * is obtained by combining (3.32), (3.37), and (3.36). It is clear 
from the form of these equations that there exists a nondecreasing function N 2 (t) such that 



\\K>i i e\\At)^N 2 (t)\\0\\(t). 
Consequently, we can write for the solution of the traction problem 

n(p, t) = w(p, t)+ J (-l)*Ef l (BEt l )*(H-Ew-Bw)(p, t). 



(3.38) 



(3.39) 



(iii) The Mixed Problem 

This problem is somewhat simpler than (ii). Thus, in the remainder of this paper we return to 
the earlier notation (2.9), (2.10). 

We seek a function u{p, t) which satisfies 



Eu+Bu=H 

0=u(b, t). 
In order to deal with the inhomogeneity in (3.40), we again assume a solution of the form 



where 



u — v-\~w 



Ev + Bv = H — Ew — Bw 



= v(b, t) 

^./ ^ T, fdVU\ t x , T r (dW . W\ , v 



0=w(b, t) 
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(3.3) 

(3.40) 
(3.41) 

(3.42) 
(3.43) 
(3.44) 
(3.45) 



An appropriate w is evidently given by 

w(p, t)=(p-b)<p(t), 

where (p satisfies 

<J>(fl,0 = [y+p(l + (~)) ^)+/ ' g(a,t,t') + m(a,t,t') (l+(^ 



a — b 



<p(t f )dt'. 



Let ft 1 be that operator which assigns to a given function #(p, t) that function A(p, which 
satisfies (3.17), (3.44), and (3.45). 

We write the general solution of (3.17) in the form 



*Y) dr+cMp+ <-f. 



h(p,t)=-J b 



Condition (3.45) implies the further restriction 

Due to (3.44), c(t) must satisfy the integral equation 



(3.46) 



l+-)F«, + 2^ 



[c(t')](t)=- z V^ 

cr ' 



JW)^ 



(0 



+ F&+jf 






0(y, *')<*y 



(t) (3-47) 



The existence of TVs (0 is clear from (3.46), (3.47), so that the solution to the mixed problem takes 
the form 

«(p,0 = w(p,0 + S (-l) n E^(BE^) n (H-Ew-Bw)(p, t). 
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